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Abstract. We give a 0(n)-time algorithm for determining whether trans¬ 
lations of a polyomino with n edges can tile the plane. The algorithm is 
also a 0(n)-time algorithm for enumerating all regular tilings, and we 
prove that at most 0{n) such tilings exist. 


1 Introduction 

A plane tiling is a partition of the plane into shapes each congruent to a fixed set 
of tiles. As the works of M. C. Escher attest, plane tilings are both artistically 
beautiful and mathematically interesting (see [2D] for a survey of both aspects). 
In the 1960s, Golomb [3] initiated the study of polyomino tiles: polygons whose 
edges are axis-aligned and unit-length. 

Building on work of Berger [2], Golomb [9j proved that no algorithm exists 
for determining whether a set of polyomino tiles has a plane tiling. Ollinger m 
proved that this remains true even for sets of at most 5 tiles. It is a long-standing 
conjecture that there exists an algorithm for deciding whether a single tile admits 
a plane tiling (see [TO:, fTT| ) 

Motivated by applications in parallel computing, Shapiro [51: studied tilings 
of polyomino tiles on a common integer lattice using translated copies of a 
polyomino. For the remainder of the paper, only these tilings are considered. 
Ollinger OH proved that no algorithm exists for determining whether sets of at 
most 11 tiles admit a tiling, while Wijshoff and van Leeuwen [55] obtained a 
polynomial-time-testable criterion for a single tile to admit a tiling. Beauquier 
and Nivat HUT] improved on the result of Wijshoff and van Leeuwen by giving 
a simpler criterion called the Beauquier-Nivat criterion. 

Informally, a tile satisfies the Beauquier-Nivat criterion if it can be sur¬ 
rounded by copies of itself (see Figure [l]). Such a surrounding must correspond 
to a regular tiling (also called isohedral ) in which all tiles share an identical 
neighborhood. Using a naive algorithm, the Beauquier-Nivat criterion can be 
applied to a polyomino with n vertices in 0(n 4 ) time. 

The 0(n 4 ) algorithm of jT] is implicit; the main achievement of [T] is a con¬ 
cise characterization of exact tiles, akin to Conway’s criterion (see usd- Gambini 
and Vuillon [8] gave an improved 0(n 2 )-time algorithm utilizing structural and 
algorithmic results on words describing boundaries of polyominoes. Around the 
same time, Brlek, Provengal, and Fedou mm also used a word-based approach 



Fig. 1 . A polyomino tile (dark gray), a surrounding of the tile (gray), and the induced 
regular tiling (white). 


to achieve 0(n )-time algorithms for two special cases: (1) the boundary con¬ 
tains no consecutive repeated sections larger than 0(y/n), and (2) testing a 
restricted version of the Beauquier-Nivat criterion (surroundable by just four 
copies). Provencal I18| further improved on the algorithm of Gambini and Vuil- 
lon for the general case, obtaining 0(n log 3 (n)) running time. In a recent survey 
of the combinatorics of Escher’s tilings, Blondin Masse, Brlek, and Labbe m 
conjecture that a 0(ro)-time algorithm exists. In this work, we confirm their 
conjecture by giving such an algorithm (Theorem [2]). 

The algorithm doubles as an algorithm for enumerating all surroundings (reg¬ 
ular tilings) of the polyomino. As part of the proof of the algorithm’s running 
time, we prove a claim of Provenqal [f§l that the number of surroundings of a 
tile with itself is 0(n) (Corollary [l]). This complements the tight bounds on a 
special class of surroundings by Blondin Masse et al. mm, and proves that 
our 0(n + fc)-time algorithm for enumerating all k surroundings (Lemma |10| is 
also a 0(n)-time algorithm. 

2 Definitions 

Here we give precise formulations of terms used throughout the paper. The 
definitions are similar to those of Beauquier and Nivat |I] and Brlek et al. 0]. 

2.1 Words 

A letter is a symbol x G £ = {u, d, 1, r}. The complement of a letter x, written 
x, is defined by the following bijection on S: u = d, r = 1, d = u, and i = r. 

A word is a sequence of letters and the length of a word W, denoted \W\, 
is the number of letters in W. For an integer i € {1,2,, \W\}, W[i\ refers to 
the ith letter of W and W[—i\ refers to the *th from the last letter of W. The 


2 


notation l k or W k denotes the word consisting of k repeats of a letter l or word 

W, respectively. 

There are several functions mapping a word IT to another word of the same 
length. The complement of W, written IT, is the word obtained by replacing each 
letter of W with its complement. The reverse of W, written IT, are the letters 
of W in reverse order. The backtrack of IT, written IT, is defined as W = W. 
Note that for any two words X and Y, AB = BA. 

2.2 Factors 

A factor of W is an occurrence of a word in W, written X A W. For integers 
1 < *, j < \W\ with i < j, W[i..j] denotes the factor of W from W[i] to W\j\, 
inclusive. A factor X starts or ends at W[i] if W[i\ is the first or last letter of 

X, respectively. 

Two factors X, Y A IT may refer the same letters of W or merely have the 
same letters in common. In the former case, X and Y are equal , written X = Y, 
while in the latter, X and Y are congruent , written X = Y. For instance, if 
W = uuulruuu then IT[1..3] = IT[6..8]. A factorization of IT is a partition of 
W into consecutive factors F\ through Fk , written W = FiF 2 ... -F*. 

2.3 Special words and factors 

A word X is a prefix or suffix of a word IT provided IT = XU or IT = UX, 
respectively. A word X is a period ofW provided \X\ < \ W\ and IT is a prefix 
of X k for some k > 1 (introduced in |T3j). Alternatively, X is a prefix of IT and 
W[i\ = W[i + \X\] for all 1 < * < \W\ -\X\. 

A factor X AW is a prefix if X starts at W[ 1], written X ^ pre IT. Similarly, 
X A W is a suffix if X ends at W[— 1], written X ^ su ff IT. A factor X AW 
that is either a prefix or suffix is an affix , written X IT. A factor X A IT 
that is not an affix is a middle , written X ^ m id IT. 

The factor X AW such that IT = UXV, \U\ = \V\, and |X| £ {1,2} is 
the center of IT. A factor X A W is a mirror, written X ^ m i r IT, provided 
IT = XUYV with Y = X and \U\ = |T|. For any X ^ m ; r IT, X refers to the 
factor Y in the definition. 

A mirror factor is admissible provided t/[l] U[— 1], T[l] V[— 1]. Observe 

that each admissible factor is the maximum-length mirror factor with its center. 
Thus any two admissible factors have distinct centers. 

2.4 Polyominoes and boundary words 

A cell is a unit square with lower-leftmost vertex (x , y) £ Z 2 and remaining 
vertices ( x + 1, y), (x, y + 1), (x + 1, y + 1). A polyomino is a simply connected 
union of cells whose boundary is a simple closed curve. 

The boundary of a polyomino consists of cell edges. The boundary word of 
a polyomino P, denoted B(P), is the circular word of letters corresponding to 
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the sequence of directions traveled along cell edges during a clockwise traversal 
of the polyomino’s boundary (see Figure [2]). 

Boundary words are circular: the last and first letters are defined to be 
consecutive. Thus for any indices i.j € Z \ {0}, W[i\ and W[i..j] are defined. 
For the boundary word W = urrdll, W[10] = W\— 9] = d and W[6..2] = lur. 



Fig. 2. A regular tiling (left) and non-regular tiling (right) of a polyomino with bound¬ 
ary word ururdrurd li luldlul. The copies in the regular tiling have a common neigh¬ 
borhood factorization ABC ABC , with A = u, B = ru, C = rdrurd. 


2.5 Tilings 

For a polyomino P, a tiling of P is an infinite set ST of translations of P, called 
copies , such that every cell in the plane is in exactly one copy. A tiling is regular 
(e.g. isohedral) provided there exist vectors o,u,v such that the set of lower- 
leftmost vertices of copies in the tiling is o + {iu + jv :i,j£ Z}. Two tilings T 
and T' are equal provided there exists a vector v such that T' = v + T. 

Copies of a tiling intersect only along boundaries, and copies with non-empty 
boundary intersection are neighbors. Lemma 3.5 of |22| implies that the inter¬ 
section between a pair of neighbors corresponds to a neighbor factor of each 
neighbor’s boundary word and these factors form a neighborhood factorization. 
Every regular tiling has a neighbor factorization common to all copies in the 
tiling. 

3 The Beauquier-Nivat Criterion 

Recall that X is the reverse complement of X. Thus X is the same path as X but 
traversed in the opposite direction. So any pair of factors X and X appearing 
on the boundary of a polyomino are translations of each other with the interior 
of the boundary on opposite sites. Beauquier and Nivat [T] gave the following 
criterion for determining whether a polyomino tile admits a tiling: 
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Definition 1. A factorization W = ABC ABC of a boundary word W is a BN 
factorization. 

Lemma 1 (Theorem 3.2 of (TV). A polyomino P has a tiling if and only if 
B{P) has a BN factorization. 

As seen in Figure [3j a BN factorization corresponds to the neighborhood 
factorization of a regular tiling. We prove this formally by reusing results from 
the proof of Lemma [T| 



Fig. 3. BN factorizations (left) and the regular tilings induced by these factorizations 
(right). For one polyomino (bottom), two of the factors are zero length. However, no 
BN factorization can have more than two length-0 factors. 


Lemma 2 (Corollary 3.2 of (TJ). Let P be a polyomino. There exists a fac¬ 
torization B(P) = F 1 F 3 F 2 F 1 F 3 F 2 if and only if there exists a tiling Sf of P with 
three copies P±, P2, P3 such that: 

— P \, P 2 , -P3 appear clockwise consecutively around a common point q. 

— Fi is the last neighbor factor of Pi whose clockwise endpoint is incident to q. 

Lemma 3. Let P be a polyomino. A factorization of B(P) is a BN factorization 
if and only if a regular tiling of P has this neighbor factorization. 

Proof. The factorization B{P) = F 1 F 3 F 2 F 1 F 3 F 2 is a generic BN factorization. 
So it suffices to prove that there exists a tiling 8 T of P satisfying the condi¬ 
tions of Lemma [2] if and only if there exists a regular tiling ^ e g with neighbor 
factorization B(P) = F 1 F 3 F 2 F 1 F 3 F 2 . 
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Tiling => neighbor factorization. Let J be a tiling and Pi,P 2 ,P 3 £ S? 
be copies as defined in the statement of Lemma [2j Let u and v be the amount 
P 2 and P 3 are translated relative to Pi, respectively. Lemma 3.2 of [TJ states 
that the copies obtained by translating Pi by u,v,v — u, —u, —v. and u — v is a 
surrounding of P\: a set of interior-disjoint copies such that every edge of Cj is 
shared by a copy. Since P 3 is a copy of P 2 translated by v — it, the neighbor factor 
of Pi incident to the__copy translated by v — u is F 2 . By similar reasoning, Pi 
has neighbor factors Pi, P 3 , and F 2 incident to the copies translated by —u, —v, 
and it— v, respectively. So Pi has neighbor factorization B(P) = FiF 3 F 2 F 3 F 3 F 2 . 
Corollary 3.1 of [T] states that for every surrounding, there exists a regular tiling 
of P containing the surrounding and thus has the neighbor factorization of Pi. 

Tiling 4= neighbor factorization. Now suppose there exists a regular 
tiling ^ eg of P with neighbor factorization FiF 3 F 2 FiF 3 F 2 . Let Pi £ be a 
copy and q be the clockwise endpoint of the factor Pi of Pi. Let P 2 , P 3 £ d^ eg 
be copies adjacent to Pi and incident to factors Pi and P 3 of Pi. Let u and v 
be the amount P 2 and P 3 are translated relative to Pi, respectively. Then q is 
the clockwise endpoint of the factor F 2 of Pi, translated by u. Also, q is the 
clockwise endpoint of the factor P 3 translated by, translated by v. So the factors 
of P 2 and P 3 whose clockwise endpoints are q are F 2 and P 3 , respectively. □ 


4 A Bound on the Number of Factorizations 

Here we prove that the number of BN factorizations of the boundary word of 
an n-omino is 0(n). This fact is used in Section [4] to improve the bound on the 
running time the algorithm from 0(n + k ) to O(n). 

Lemma 4. Let W be a boundary word with a factor X. Let P, S P m i r W such 
that P ^ pre X, S ^suB X, and P ^ S. Then X has a period of length 2|Aj — 

Proof. Since P and S are mirror, there exists X' <W with \X'\ = |Aj, P ^ pre 
A'', and S ^ S uff X'. Observe that X has a period of length r > 1 if and only 
if X[i\ = X[i + r] for all 1 < i < \X\ — r. Let 1 < i < |P| + |S| — |Aj. Then 
1 < \P\ + 1 - i < | A'| and 1 < \P\ + 1 + |5| - \X'\ -i<\S\. So: 

AW = P[i] 

= P[|P| + l-<] 

= ^[|P| + 1-*] 

= f[|P| + l + |5|-|A / |-i] 

= S[\S\ + l-(i + \X / \-\P\)] 

= S{ l + \X'\-\P\] 

= X[i + \X'\ — |P| + (|A| — IPI)] 

= X[i + 2\X\ — (|P| + |Sj)] 
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Since P ^ S, 2\X\ — (|P| + \S\) > 2\X\ — {2\X\ — 1) = 1. So X has a period of 
length 2|A| - (|P| + |5|). □ 


Lemma 5. Let W be a boundary word with X A W. Let P, S ^ m ir W such 
that P ^ pre X, S S su ff X, and P ^ S. Any factor Y ^ m id X with |Y| > 

2|A| — (|P| + |5|) is not an admissible factor ofW. 

Proof. By Lemma [dj X has a period of length r = 2\X\ — (|P| + IS 1 !). Let 
Y ^ mid X and \Y\ > r. 

Let X' < W with \X '| = |A| and the center of X' exactly \W\/2 letters 
from the center of X. Then P A pre X', S ^ su ff X' , and Y ^ m id X'. Again by 
Lemma |4j X' has a period of length r. 

Let U,V A W such that W = YUYV. Since Y is a middle factor of X , 
the letter C/[l] is in X. Since X has a period of length r and |Y| > r, Z7[l] = 

Y[|Y| + 1 — r] = Y[r]. Since Y is a middle factor of X' and X' has a period of 
length r, U[— 1] = Y[r\. So U[ 1] = U[— 1] and Y is not admissible. □ 

Lemma 6. Let W be a boundary word. There exists a set LX o/0(l) factors of 
W such that every F ^ a dm W with |F| > |W |/6 is an affix factor of an element 
of-?. 

Proof. A special case on three factors. Let Pi,P 2 ,P 3 ^ a dm W with |Pi|, |P 2 |, |P 3 | > 
|W |/6 and centers contained in a factor of W with length at most |W|/14. Let 
A' A W be the shortest factor such that Pi,P 2 ,P 3 A A, and so Pi ^ pre X and 
Pj ^suff X for some i, j £ {1, 2 ,3}. We prove that if i ^ j. then Pi, P 2 , P 3 ^ a ff A. 

Without loss of generality, suppose i = 1, j = 2 and so P 3 ^ m id A. By 
Lemma [5j since P 3 ^ adm W, |P 3 | < 2 |A| - (|Pi| + |P 2 |) < |P 3 | + \W\/7 + |P 2 | - 
(|Pi| + |P 2 |) = |W|/7 < |W|/ 6 , a contradiction. So P 3 ;< a ff X. 

All nearby factors. Consider a set J? = {Pl, P 2 ,..., F m } of at least three 
admissible factors of W of length at least |W |/6 such that the centers of the 
factors are contained in a common factor of W of length |W|/14. We will prove 
that every element of J? is an affix factor of one of two factors of W. 

Let G A W be the shortest factor such that Pj A G for every Pj £ jP It is 
either the case that there exist distinct Fi,F r £ J* with Pj ^ pre G, F r ^ su ff G, 
or that G £ JL and every Pj £ J? besides G has Pj ^ m id G. 

In the first case, Pj ;< a ff G for any j / !, r by the previous claim regarding 
three factors. Also Pj,P r ^ a s G. So every factor in J? is an affix factor of G. 

In the second case, let G' A G be the shortest factor with the same center as 
G such that every factor in excluding G is a factor of G'. Clearly G' ^ m ir W 
and G' ^ a dm W. Without loss of generality, there exists F p £ ^ such that 
F p ^ pre G'. Since F p ^ adm W and G' Adm W, F p ± G'. 

Applying Lemma [5] with X = G'. P = F p , S = G', every middle factor of G' 
in X has length at most 2\G'\ — (|G'| + |F p |) < \G'\ — |P p | < |W|/7 < |W|/ 6 . So 
every factor of G' in J? is an affix factor of G'. Thus every factor in J? is either 
G or an affix factor of G'. 

All factors. Partition W into 15 factors p, p,..., p 5 each of length at most 
| IP |/14. Let be the set of admissible factors with centers containing letters in 
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Ii. Then by the previous claim regarding more than three factors, there exists a 
set (G and possibly G') such that every element of J'i is an affix factor of an 
element of and |^| < 2. So every F ^ a d m W with |F| > |W|/6 is an affix 
factor of an element of AA = (J^ and \&\ <2-15. □ 

Theorem 1. A boundary word W has 0(|W|) BN factorizations. 

Proof. Consider the choices for the three factors A, B , C of BN factorization 
W = ABC ABC. In any factorization, some factor has size at least |W|/6. By 
Lemma [hj there exists a 0(l)-sized set of factors & such that any factor with 
length at least |W|/6 is an affix factor of an element of &. Without loss of 
generality, either \A\ > |W|/6 and A is a prefix of a factor in & or \C\ > |W|/6 
and C is a suffix of a factor in &. 

Let H = ABC be the factor formed by consecutive factors A, B, C of a BN 
factorization. Then since \H\ = |W|/2 and shares either the first or last letter 
with a factor in jF, there are 0(1) total factors H. For a fixed H, choosing the 
center of B determines B (since B is admissible) and thus A and C. So there 
are at most 2(|W|/2) factorizations for a fixed factor H. □ 

Since Lemma [3] proves that factorizations and tilings are equivalent, the pre¬ 
vious theorem implies a linear upper bound on the number of regular tilings of 
a polyomino: 

Corollary 1. An n-omino has 0(n) regular tilings. 

As pointed out by Provengal [H], it is easy to construct polyominoes with 
Q(n) such tilings. For instance, the polyomino with boundary word W = ur'dF 
with i > 1 has |W|/2 — 1 regular tilings. 


5 An Algorithm for Enumerating Factorizations 

The bulk of this section describes a 0(|fF|)-time algorithm for enumerating 
the factorizations of a polyomino boundary word W. The algorithm combines 
algorithmic ideas of Brlek, Provengal, and Fedou [2] and a structural result based 
on a well-known lemma of Galil and Seirferas [5]. 

Lemma 7 (Corollary 5 of |4j)- Every factor of a BN factorization is admis¬ 
sible. 

Lemma[8]is a variation of Lemma C4 of Galil and Seirferas j5|. We reproduce 
their proof with minor modifications. 

Lemma 8. Let A and B be two words of the same length. Moreover, let A = 
X 1 X 2 = YiT 2 = E\’Z 2 and B = X Q X 2 = Y& = Z X Z Q with \X x \ < |W| < \Z X \. 


A 


x 2 

Y 1 

Y 2 


V 


Zi 

Z2 


B 

Xq 

x 2 

Y 1 

y 2 


V 


Zr 


Fig. 4. The words used in the proof of Lemma [8] 


Proof. Let V be the word such that \\V = Z\ (see Figure [4]). 

Claim (1): V is a period of Z-\ . Since Y{V = Z\ , then Z\ = Y{V = VY\ 
is a prefix of B. So_Yi is_a prefix of Z\ = VYi and thus V is a period of Y\. So 
V is a period of VY 1 = Z\. ^ 

Claim (2): V is a prefix of X 2 . Since V is a prefix of Y 2 , V is a suffix of 
Y 2 . So V is a suffix of X 2 and V is a prefix of X 2 . 

Claim (3): X\V is a prefix of Z\. Since V is a prefix of X 2 , XiV is a 
prefix of Y{V. Since |-XiV| < \Y 1 V\ = \Z 1 \, X{V is also a prefix of Z\. 

Claim (4): V is a period of X\. By claim (1), V is a period of Z\ 1 so Z\ 
has a period of length |1/| = \V\. By claim (3), X\V is a prefix of Z\ and so also 
has a period of length \V\. Then X\V = VX\ has a period of length \V\, namely 
V. So V is also a period of X-\. 

Finally, combining claims (1) and (4), since V is a period of both Xq and 
Xi, Xq = X±. By symmetry, the same proof also implies Zq = Z^. □ 

Lemma 9 (Theorem 9.1.1 of |12] h Two non-circular words X, Y can be 
preprocessed in 0(\X | + |F|) time to support the following queries in 0(l)-time: 
what is the longest common factor of X and Y starting at X[i\ and Y[j ] ? 

Lemma 10. Let W be a polyomino boundary word. Then the BN factorizations 
ofW can be enumerated in 0(|W|) time. 

Proof. Lemma [7] states that BN factorizations consist entirely of admissible fac¬ 
tors. The algorithm first computes all admissible factors, then searches for fac¬ 
torizations consisting of them. 

Computing admissible factors. Lemma [T] implies that there are at most 
2\W\ admissible factors, since admissible factor has a distinct center. For each 
center W[L.i] or W[i..i + 1], the admissible factor with this center is LR, where 
R is the longest common factor of W starting at W[i + 1] and W starting at 
W r [|l / F|/2 — (i + 1)]. Similarly, L is the longest common factor of W starting at 
VF[|VF|/2 — i] and W starting at W[i\. Preprocess WW, WW, WW, and WW 
using Lemma 9] so that each longest common factor can be computed in 0(1) 
time. If \L\ ^ 1?|, then X is not admissible and is discarded. Since 0(1) time is 
spent for each of 2\W\ admissible factors, this step takes 0(|W|) total time. 
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Enumerating factorizations. Let W = AYAZ with A an admissible factor 
and \Y\ = \Z\. Let Bi, B 2 , ■. ■, Bi be the admissible prefix factors of Y, with 
|Bi| < \B 2 \ < ••• < \Bi\. Similarly, let C\,... ,C m be the suffix factors with 
|Ci| < ••• < \C m \. Lemma [8] implies that for fixed A 1 there exist intervals 
[b, l ], [c, to] such that the BN factorizations ABiCjABiCj are exactly those with 
i £ [b , l] or j £ [c, to] . 

First, construct a length-sorted list of the admissible factors starting at each 
W\k\ in 0(|W|) time using counting sort. Do the same for all factors ending at 
each W[k}. 

Next, use a two-finger scan to find, for each factor A that ends at W[k\, the 
longest factor Bi starting at W[k + 1] such that |A| + \Bi\ < |W|/2. Then check 
whether Cj, the factor following Bi such that \ABiCj\ = |W|/2, is admissible and 
report the factorization ABiCjABiCj if so. Checking whether Cj is admissible 
takes 0(1) time using an array mapping each center to the unique admissible 
factor with this center. 

Additional BN factorizations containing A are enumerated by checking fac¬ 
tors Bi with i = l — 1, l — 2,... for an admissible following factor Cj. Either Cj 
is admissible and the factorization is reported, or i = b — l and the iteration 
stops. 

Finally, use a similar two-finger scan to find, for each factor A that starts 
at W[k\, the longest factor C m that ends at W[k + |W|/2 — 1] such that |A| + 
\C m \ < |IF|/2, check whether Bi preceeding C m such that \ABiC m \ = |W|/2 
is admissible, and report the possible BN factorization. Then check and report 
similar factorizations with Cj for j = m — 1, to — 2,... until j = c — 1. 

In total, the two-finger scans take 0(|W|) time plus 0(1) time to report 
each factorization. Reporting duplicate factorizations can be avoided by only 
reporting a factorization if A[l] appears before B[ 1], C[l], A[l], B[ 1], and 0[1] 
in W. Then by Theorem]!] reporting factorizations also takes 0(|W|) time. □ 

Combining this algorithm with Lemmas[l]and[3]yields the desired algorithmic 
result: 

Theorem 2. Let P be a polyomino with n edges. In 0(n) time, it can be deter¬ 
mined if P admits a tiling and the regular tilings of P can be enumerated. 
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